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General quantum-mechanical description of relativistic particles and nuclei with spin 1/2 
channeled in bent crystals is performed with the use of the cylindrical coordinate system. The 
previously derived Dirac equation in this system is added by terms characterizing anomalous 
magnetic and electric dipole moments. A transformation to the Foldy-Wouthuysen represen¬ 
tation, a derivation of the quantum-mechanical equations of motion for particles and their 
spins, and a determination of classical limit of these equations are fulfilled in the general 
case. A physical nature of main peculiarities of description of particles and nuclei in the 
cylindrical coordinate system is ascertained. 

1 Introduction 

The strict quantum-mechanical description of relativistic spin-1/2 particles and nuclei chan¬ 
neled in bent crystals is an important problem [T], Taking spin effects into account is very 
important in such a description. It is known that, during planar channeling in bent crystals, 
the particle and nucleus spins are rotated by a rather large angle. This effect was first found 
in the papers of Baryshevsky M, where he also proposed its use to determine the mag¬ 
netic moments of shortliving particles. The simple dependence between the rotation angles 
of particles and their spins in bent crystals was found by Lyuboshits jlj. In such crystals, the 
centrifugal force acting on particles or nuclei moving along bent trajectories is compensated 
by the Coulomb force, which leads to the appearance of a rather strong electric field rotating 
the spin. The effect of spin rotation was observed experimentally in [5, 6]. 

Although particle and nucleus channeling in many cases can be adequately described 
by methods of classical theory, a thorough quantum-mechanical analysis of the problem 
is also necessary. So, the discreteness of the energy spectrum is very often important for 
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relativistic positrons (electrons). To determine a given spectrum correctly, it is necessary 
to adequately take into account spin effects for relativistic particles. As is well known, in 
the Dirac equation, as in the Dirac-Pauli equation, which takes the anomalous magnetic 
moment (AMM) into account, Dirac matrices determine the spin projections on the axes of 
the Cartesian system of coordinates. It is convenient to use such a system of coordinates 
only for channeling in unbent crystals, and the choice of the cylindrical system of coordinates 
is natural for bent crystals (if the radius of curvature is approximately constant). 

The author of [7] presented the quantum-mechanical description of spin 1/2 particles 
(nuclei) for planar channeling in unbent and bent crystals. Special attention was focused 
on determination of the spin dynamics. In this paper, the author solved the Dirac equation 
(supplemented with terms describing the AMM) in the Foldy-Wouthhuysen (FW) represen¬ 
tation, constructed the operator equation of spin motion, and calculated the average value 
of its precession frequency. The results obtained using the quantum-mechanical description 
agreed entirely with the corresponding classical results. The author of [7] used the Carte¬ 
sian system of coordinates rather than the cylindrical one and took the presence of crystal 
bending into account formally by including an additional potential energy determining the 
correction for the centrifugal force into the Hamilton operator in the FW representation. 

Naturally, such an approach is not rigorous, although it leads to reasonable results. In 
the present paper, the Dirac equation in the cylindrical system of coordinates derived in [8j 
is used as the initial one. We supplement it with terms describing the AMM and the electric 
dipole moment (EDM) and perform transformation into the FW representation by a method 
developed in [9]. We use the obtained Hamilton operator in this representation to derive 
general equations describing the motion of particles and nuclei and the spin evolution. 

We let the respective Greek and Roman letters = 0,1, 2, 3 and i, j, k,... = 

1, 2, 3 denote the world and space indices in four-dimensional spacetime. Using the apparatus 
of the theory of general relativity, tetrad indices are denoted by the initial letters of the 
Roman alphabet a, 6, c,... = 0,1,2,3. The time and space tetrad indices are singled out 

by hats. The signature has the form (H-). We here use the system of units fi = 

1, c = 1. In some cases, to make the presentation clearer, we include Plancks constant in 
the corresponding formulas. The notations determine the commutators 

and anticommutators, respectively. 
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2 Dirac-Pauli equations in the cylindrical system of 


coordinates 


The standard DiracPauli equation (in the Cartesian system of coordinates) has the form 


7 m tt m -m + 


T = 0, 


( 1 ) 


where and a = i (y^y^ — 7 I/ 7 /i )/2 are the Dirac matrices, F fJll = (E,B) is the electro¬ 
magnetic field tensor, p' is the AMM, and 717 = iD p = ih(d/dx M ) — eA M . Here .E, B and 
are the electric strength, the magnetic induction, and the four-potential of the electro¬ 
magnetic field. In [lOj, this equation was supplemented with a term describing the EDM 


d: 
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m + ^■cr fJ ’ l 'F flv 
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fllS 


T = 0, 


( 2 ) 


where G^ = (-B, E ) is the tensor dual to F^. 

In principle, it is possible to pass to cylindrical coordinates, not changing the definition 
of the Dirac matrices and using the relations 


7 P = 7 ■ e p = 7 1 cos (j) + "f y sin 0 , 7 ^ = 7 • e<j, = —y 1 sin 0 + 7 ^ cos 0 . 

Naturally, such a way is not convenient, and the subsequent transformations are accompanied 
by very cumbersome calculations. Of course, all required calculations can be carried out in 
the Cartesian system of coordinates. However, to specify an external field, if the symmetry 
of the problem is taken into account, it is more convenient to use the cylindrical coordinates 
specifically. A convenient form of the Dirac equation in the cylindrical system of coordinates 
was found in [ 8 ], the authors of which used the fundamental Pauli theorem HU determining 
the relationship between different sets of Dirac matrices satisfying the required commutation 
and anticommutation relations. The Dirac equation derived in [ 8 j (for p = 0) has a very 
simple form and formally coincides with the initial equation: 


(7^77 — m)T = 0. 


(3) 


The matrices y^ are ordinary Dirac matrices. However, here, the indices 1, 2, and 3 corre¬ 
spond to the cylindrical coordinates p, 0 , and z, and 


(tTi, 7T 2 , 7T 3 ) = 


eA p , 

pocp 


d N 
eA#, ih— - eA z 


(4) 
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The contravariant vector potential A , whose components are —A, is usually used. 

It follows from Eq. (J3J) that the operator i q contains the nabla operator in the cylindri¬ 
cal system of coordinates. The Dirac equation in the spherical system of coordinates has 
similar properties [ 8 ]. The result obtained in | 8 ] is completely natural and shows that the 
transformation into cylindrical and spherical coordinates which retains the form of the Dirac 
matrices y^does not violate the covariance of the Dirac equation. 

We use this fact to include terms proportional to the AMM and the EDM into the 
equation. Such inclusion will be substantiated additionally in the next section by comparison 
with the results obtained within the framework of the theory of general relativity. 

The covariance of Eq. (EJ) is not violated if its generalization is written in a form (|2]h 
where the matrices have the usual form and the indices /j and v correspond to the 
cylindrical system of coordinates. The components of the tensors F and G IIU are also 
determined in this system of coordinates. 

It is convenient to multiply the obtained equation by the matrix 7 0 and represent it into 
the Hamiltonian form 

= "HT, V. = ol ■ tv + /3m + e<f> + //(—II • B + ij ■ E) — d(H ■ E + zy • B ), 

7r = — ihV — eA, 

where $ = A 0 is the scalar potential and tv = — (tt p , 7 r 2 ) is the kinetic momentum operator 
in the cylindrical system of coordinates. Unlike the Cartesian system of coordinates, the 
different components of the operators tv and V do not commute with each other in the 
general case. 

Equation (J5j) is initial for the next transition to the FW representation. Nevertheless, the 
obtained result requires additional substantiation because the terms describing the AMM 
and EDM were included in ([3]) obtained as a result of the transformation, rather than in 
the standard Dirac equation. Such substantiation can be obtained within the framework of 
gravitation theory. Using gravitation theory makes it possible to relatively simply understand 
the physical meaning of the difference between the dynamics of particles in the Cartesian 
and cylindrical systems of coordinates. 
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3 Deriving the Dirac-Pauli Hamiltonian for particles 
in the cylindrical system of coordinates using the 
methods of gravity 


The correctness of the above generalization of Eq. (j3J), which makes it possible to take into 
account the possible presence of the AMM and the EDM of the particle, can be confirmed 
by analyzing the covariant Dirac equation in gravitation theory. This equation describes the 
electromagnetic interaction of the Dirac particle in Riemannian spacetime and has the form 

(ifi^ a D a — m)ijj = 0, a = 0,1,2,3. (6) 

In Eq. (JH]) we have 

Da = < (A + jA^ + (7) 

where r abc = —Tb ac are the Ricci rotation coefficients [14]) having the following form: 

r a fe c = ^ ( C abc - C bca - C cab ), C abc = e^ebid^ - d u e cll ). (8) 


Here are tetrad coefficients determining the tetrad components of the covariant derivative 
( D a = e^D^) and other four-vectors. 

General equation ([6]) can be used for cylindrical and other curvilinear coordinates. The 
transition from x, y , z to p, 0, 2 in the expression for the squared infinitesimal interval ds 2 
gives the following form for the metric tensor: 


9 iw 


' 1 0 0 0 N 

0-10 0 
0 0 -p 2 0 

v 0 0 0 -1 ) 


(9) 


Metric tensor ([9]) describes a flat spacetime. 

In this case, the nonzero Ricci rotation coefficients are 


1^212 — — r i2 2 — -• (10) 

p 

The covariant Dirac equation obtained by substituting (TTOl) into ([6]) and (J7]) does not coincide 
with Eq. ([3]). However, we note that to find the Hermitian Hamiltonian, it is necessary to 
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perform the nonunitary transformation of the wave function of the initial covariant Dirac 
equation © which has the form [15] 

D 1 / 2 nk — 


v = ^ = y/pnl), 


( 11 ) 


where g is the determinant of the metric tensor. To obtain the Hermitian Hamiltonian from 
Eq. (|3j), such a transformation is not required. 

The natural generalization of Eq. (J6]) making it possible to describe the AMM and the 
EDM and retaining the equation covariance has the form 




T = 0. 


( 12 ) 


Here, F ab = e£e" b F pu: G ab = e^e b G pu , and o ab = i{^ a ^ b - yV)/ 2 - 

In the general case, Eq. ([ 6 ]) was transformed into the Hamilton form in [[15], and the 
authors of [16] transformed the obtained Hamiltonian into the FW representation also in the 
general case. For metric (15]) . using the method for finding the Hamilton form of Eq. (U5D 
proposed in [15] leads to an expression for the Hermitian Hamiltonian which coincides with 
(l5|) . Thus, the correctness of Eq. ([5]) is completely confirmed strictly by the methods of 
quantum mechanics of Dirac particles in gravitational fields. 

Gravitation theory also enables us to determine the distinguishing features of particle 
dynamics in the cylindrical system of coordinates. The particle velocity in this system 
has the form v = v p e p + + v z e z . Naturally, its constancy (v = const ) means that 

a force, keeping the particle in a circular orbit, acts on it. This force can be determined 
using the formalism of gravitoelectromagnetic fields which makes it possible to describe 
relativistic particles in arbitrarily strong gravitational fields. This formalism was proposed 
by Pomeranskii and Khriplovich mi and underwent further development in [T5lfl 6 l fl 8 l 19] . 
It is convenient to introduce gravitoelectromagnetic fields to determine the equations of 
motion described in local Lorentz (tetrad) reference systems: 


ds 

dt 


du 


= n x s, f 2 = 




u 


n-v»\ e+ 


uxB" 


w u 


-B + 

difi 

dt 


u x S 

vP + 1, 


S ■ u 


n u 


(13) 


(14) 


These equations are analogous to equations for a particle with a Dirac magnetic moment 
(g = 2) in an electromagnetic held. The spin s, the gravitoelectric held £, and the gravito- 
magnetic held B are determined precisely in these local Lorentz reference systems; i.e., their 
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components are of the tetrad type. The general equation for the fields has the form 

£■ = r 0 ic u c , £ 7 ■= — -e-ikiT k i c u c , (15) 

where e^-i is the antisymmetric tensor with spatial components. 

The tetrad and world directions coincide in terms of the problem under consideration, 
and the fields are determined by the expressions 

e = 0, B-=B 2 = 0 , B?=— . (16) 

Because et. = 1/p, we have u ^ = vP / p. 

Equations (THj) and (fl6|) show that the force determined by the gravitomagnetic field 
acting in the cylindrical system of coordinates is an analogue of the Lorentz force. Its 
appearance is a consequence of the fact that, if the azimuthal angle of the particle changes 
by d(f>, the horizontal axes of the cylindrical and Cartesian systems of coordinates rotate by 
the same angle with respect to each other; i.e., the cylindrical system of coordinates rotates 
with an instantaneous angular velocity —d(j)/dt = —vp/p with respect to the Cartesian one 

m- 

For the detailed quantum-mechanical description of spin 1/2 particles and nuclei, it is 
convenient to use the FW representation, the transition to which is performed in the following 
section. 


4 Dirac-Pauli Hamiltonian in the Foldy-Wouthuysen 
representation for particles in the cylindrical system 
of coordinates 

The method of FW transformation for relativistic particles in arbitrary external fields de¬ 
veloped in [SET, l22| 123]. is based on the following representation of the initial Hamiltonian: 

U D = /3m + £ + O, P£ = SP, PO = -OP, (17) 

where P = 7 0 . I 11 this equation, the initial Hamiltonian is divided into even (£) and odd 
(■ O ) terms, which are diagonal and off-diagonal in two spinors, respectively. I 11 Eq. (fTTlh we 
have 

£ = e<h — p 1 II • B — dll ■ E, O = a ■ 7r + ■ E — id"y ■ B. (18) 
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The FW transformation proved itself to be the best method for finding the quasiclassical 
approximation and the classical limit of relativistic quantum mechanics in the case of the 
single-particle approach. In this representation, the Hamiltonian is diagonal in two spinors 
(block-diagonal), the probability interpretation of the wave function is restored, and the 
operators have the same form as in nonrelativistic quantum mechanics. The transition to 
the FW representation is widely used for all fundamental interactions. 

If the terms that are bilinear in external fields E and B are disregarded, the result of 
the transformation can be written in the form 

H FW = Pe + S-\{ -t [0,[0,J]]|, t = VF? + O 2 . (19) 

8 (e(e + m) ) 

This Hamiltonian in the FW representation contains exactly determined terms of the first 
and second orders with respect to h/S 0 , where S 0 is the value of the action dimensionality 
[23] . Terms of second and higher order with respect to h/So are also determined exactly in 
the case where they appear as a result of calculating the given Hamiltonian. In particular, 
this is related to the Darwin (contact) interaction. 

For the problem under consideration, the commutator of the operators 7r* and tt j is equal 
to the sum of the two terms: 


[vh, 7 Tj] = -h 2 [Vi, Vj] - ieh(yiAj - VjAi) = - h 2 [\/i , Vj] + iehe ijk B k . 

The first term is equal to zero if Cartesian coordinates are used. For cylindrical coordinates, 
the commutator 

has a nonzero value. 

As a result, the operator O 2 is defined by the following exact expression: 

O 2 = tt 2 - ehS • B - hzM + 0 (S • [tt x G\ - S • [G x tt] - W • G) + G 2 , (20) 

P 

where G = n'E — dB. 

When calculating the general expression for the Hamiltonian in the FW representation, it 
is possible to disregard terms of second and higher degrees with respect to external fields and 
terms of third and higher degrees with respect to the Planck constant. With this accuracy, 
the Hamiltonian is determined by the equations 

nj _ -i/(0) I nj(MDM) ^j(EDM) /<->-, \ 

ttFW — ft-FW ' FW ' 'V FW ’ 




( 22 ) 
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p 0 m , 
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e + m 


'Z-lnxE] - H-\Ex ttI - ftV-B 

ns 


+/4 


4 1 e'(e' + m) ’ . 


'S-7r)(S-7r) + (S-7r)(7r-S) + 27rft(7r-j + j ■ 7r 


(23) 


M 


[EDM) _ 


FW 


= -dH-E + 

_ d 
~4 


d 


4 1 e'(e' + m) ’ . 


^•7r)(n-7r) + (n-7r)(7r-^) 


'E-lnxB] — HABxir] 


(24) 


where po = eh/{2m) is the Dirac magnetic moment and j is the density of the external 
current satisfying the Maxwell equation 


3 = 


47r 


V x B 


dE \ 

~dt J ‘ 


The quantities ft^P M ^ and define the contributions of the magnetic and electric 

dipole moments (the MDM and the EDM), respectively. In accordance with the Maxwell 
equations 

dB 

V ■ -B = 0, V x E = ——, 

the EDM does not contribute to contact interactions with external charges and currents. 

Equations give a general solution of the quantum-mechanical description of 

the electromagnetic interaction of the Dirac particle in the cylindrical system of coordinates. 


Comparison with the results in [10] shows that the operators and ft 

the same form as in the Cartesian system of coordinates, and a new spin-dependent term 


I [EDM) 
FW 


have 


appears in ft 


(o) 

FW- 


5 General equations of particle dynamics 


Deriving quantum-mechanical equations of particle and nucleus dynamics in the cylindrical 
system of coordinates has important specific features related to the noncommutation of the 
operators V p and V^. The form of the dynamic equations is universal: 

dA 


dn i , cbr i 

~^ = ^fw,n] + — = -[H FW ,7r] 

c/n i. , „ 


at ’ 


(25) 

(26) 


9 












where Q is the operator of the angular velocity of the spin precession. However, the above 
noncommutation leads to the appearance of new terms in dynamic equations compared to 
the corresponding equations in the Cartesian system of coordinates. 

As is known, the total force acting on the particle with spin has a spin-dependent compo¬ 
nent. This component is usually called the Stern-Gerlach force, and, for a moving particle, 
it depends on the magnetic and electric fields. Although the Stern-Gerlach force leads to 
the important effect of division of the beam into two beams with different polarizations, as a 
rule, it is small for charged particles as compared with the Lorentz force e(E + v x B). The 
author of [7j showed that, in the case of planar channeling, the influence of spin-dependent 
terms in the Hamiltonian in the FW representation on the particle trajectory is small and 
was not observed experimentally. Therefore, we do not take the spin-dependent force into 
account when determining the dynamics of the kinetic momentum of a particle. 

In the approximation used, the equation of motion for kinetic momentum (|25j) takes the 
form 



The first two terms in the obtained equation give the operator expression for the Lorentz 
force, and the third term describes the additional force acting in the cylindrical system of 
coordinates. It can be written in a more compact vector form. We introduce the operator 
O = Oe z = ( Pcf,/p)e z , characterizing the orbital motion of the particle about the z axis. In 
this case, 

(ttxO-Oxtt)}. (28) 

The author of [22j showed that, if the conditions of the quasiclassical approximation 
(the de Broglie wavelength is less than the characteristic size of the region of external held 
nonuniformity) are satisfied, using the FW representation it becomes possible to reduce de¬ 
termination of the classical limit of equations of relativistic quantum mechanics to replacing 
operators in the Hamiltonian and the quantum-mechanical equations of motion with the cor¬ 
responding classical quantities. In this case, it is possible to disregard the noncommutation 
of the operators in the quantum-mechanical expressions. In the case under consideration, 
O — y e'uj = mycu, where u; = c oe z = {y^/p)e z is the instantaneous angular velocity of the 
orbital particle motion about the z axis and 7 is the Lorentz factor. The classical limit of 
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the expression for the additional force has the form 


T —y —uj x 7 r. 


( 29 ) 


It is easy to see that, for particle motion along a circular arc, this equation reproduces the 


centrifugal force TT^v^/p, and the force T changes the instantaneous angular velocity of the 


orbital particle motion by the quantity — u>. 

Expression (129|) is completely analogous to the corresponding, expression for the force 


acting in a rotating reference system (Eq. (3.37) in Ref. The analogy between 


the cylindrical system of coordinates and the rotating reference system was demonstrated 
previously in [20] within the framework of the classical approach. 

However, it is important that the rotating reference system can be used only within the 
framework of the one-particle description, while the cylindrical system of coordinates can 
also be used to describe a beam of particles or nuclei with various energies. The presence of 
such a possibility is an important advantage of the cylindrical system of coordinates. 

The force T is similar to forces appearing in noninertial reference systems. However, the 
latter are real forces, which, in particular, can be measured by a dynamometer, while the force 
T is fictitious. Its presence does not affect the dynamometer readings, although it affects the 
particle and nucleus motion in the cylindrical system of coordinates. An interesting effect of 
the mutual influence of particle motion in the directions e p and follows from Eqs. (I27|) 
and (j29l) . In particular, oscillator motion in the radial direction causes the appearance of an 
oscillating term in the expression for p$ and vice versa. This effect has kinematic nature. It 
is due to the rotation of the axes e p and (about the Cartesian system of coordinates) as 
the azimuth (j) changes, and it disappears at p^ = 0 as follows from the definition of O. 

Spin motion in the cylindrical system of coordinates has a simpler character than the 
dynamics of the kinetic momentum. The angular velocity operator of the spin precession, 
which is easily determined using Eq. ([26]) . has the form 


n = n(°) + q ( mdm ) + EDM \ 


(30) 


where 



(31) 


1 In Ref. TB], as in Eqs. (1271) and (OBI) , the force was defined as the time derivative of the covariant 


momentum operator. Coriolis factor 2 appears in the corresponding equations for the acceleration and the 
spatial component of the contravariant four-velocity m mi¬ 
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q(edm) = _ / 3 ld E ± 

h 2 h 
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>K{(T + 4 *} 


e'(e' + m) ’ 


' E ■ 77)77 + 77(77 • E ] 


— < —, (77 x B — B x 77' 

27z'-' ,v 


( 32 ) 


(33) 


For a particle with a positive total energy in the FW representation, the lower spinor is zero. 
Therefore, passing to the classical limit eliminates the (3 matrices (and also the anticommu¬ 
tators). 

It is easy to see that Q , (0 ^ —>■ — uj in the classical limit. Thus, the passage from the Carte¬ 
sian to the cylindrical coordinates changes the instantaneous angular velocities of orbital 
particle motion and the precession of its spin by the same quantity, so that the difference 
between them remains the same. The forms of the quantities £l (AfDM ' 1 and Q, {E,)hl> in the 
Cartesian [10] and cylindrical systems of coordinates coincide. 


6 Motion of particle and nucleus spins at planar chan¬ 
neling in bent crystals 

As was mentioned above, the quantum-mechanical description of spin 1/2 particles and 
nuclei during planar channeling in unbent and bent crystals was presented previously in [7], 
the author of which used the Cartesian system of coordinates, and the presence of crystal 
bending was taken into account by including an additional potential energy into the Hamilton 
operator in the FW representation. In our notation, this energy is 

W = (34) 

li 

where R is the crystal bending radius and x = 0 corresponds to the middle of the distance 
between the crystal planes. During planar channeling in bent crystals, the interplanar dis¬ 
tance d p is negligibly small as compared with the crystal bending radius. Therefore, the 
approach used in [7] leads to the obtained correct results. The author of this paper gave 
detailed quantum-mechanical description of the effects occurring during particle and nucleus 
channeling in unbent and bent crystals. For this reason, we restrict ourselves to deriving 
well-known formulas M, describing the spin motion during planar channeling in bent 
crystals, starting from the general equations obtained above. 
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At planar channeling, the held of planes is characterized by the even potential $(x) = 
$(—x). For nuclei that move in the channeling mode and have positive charges, this held 
can be approximated by the harmonic potential 


T , ax 2 817o 

= “= IF' 


(35) 


where Uq is the maximum value of the potential and d p is the distance between the crystal 
planes. We assume that the crystal is bent so that the bending plane is perpendicular to the 
crystal planes and the curvature radius is R. In this case, the plane potential (135]) becomes 

a(p — R) 2 


*(p) = 


(36) 


We do not consider the magnetic held and neglect effects caused by possible electric dipole 
moments. 

At channeling, the particle oscillates with respect to the equilibrium trajectory which is 
a circular arc. Therefore, the average force acting on the particle in the cylindrical system 
of coordinates is zero. As follows from Eqs. (127|) and (1291) . in this case, in the classical limit, 
we have 


e < E >=< w x 7r > . 


Another form of this equation is 


e < E p >= 


— m\u 



- 1 . 


(37) 


The fact that the kinetic energy of transverse motion is significantly smaller than that of 
longitudinal motion during channeling is taken into account in Eq. (I37lh 

Substitution into Eqs. (130]) and (132|) gives the angular velocity of the spin precession 


n 


l \g-2 
7 l 2 



<JL>. 


(38) 


When passing to the Cartesian system of coordinates, the quantity —a; is added to the an¬ 
gular precession velocity, and we obtain the Lyuboshits formula [I] determining the relation 
between the angular velocities of the spin and momentum rotations: 


j^(Car) 


7 ~ 1 

7 


9 ~2 


(7 + 1) + 1 


U). 


(39) 


The formula determining the effect of particle and nucleus spin rotation during planar 
channeling in bent crystals was first derived in a paper by V.G. Baryshevsky [2j. The fast 
increase in the ratio Vt/uj with increasing Lorentz factor makes it possible to measure the 
magnetic moments of relativistic particles with short lifetimes [2], 
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7 Discussion and summary 


It is natural to choose the cylindrical system of coordinates in the quantum-mechanical de¬ 
scription of relativistic spin 1/2 particles and nuclei channeled in bent crystals. However, the 
problem of determining projections of the spin operator on the radial and azimuthal direc¬ 
tions appears in this case. Using ordinary Dirac matrices for the given projections is its best 
solution, found in [8j. In this paper, the Dirac equation in cylindrical coordinates obtained 
in [8] is supplemented with terms describing the AMM and EDM (with strict substantiation 
of this procedure using methods of quantum mechanics of Dirac particles in gravitational 
fields). Although the form of the obtained equation does not differ from that of the cor¬ 
responding equation in Cartesian coordinates, the difference between them is manifested 
after transformation into the FW representation. In this representation, the Hamiltonian in 
cylindrical coordinates differs from the corresponding Hamiltonian in Cartesian coordinates 
by the presence of an additional spin term. Although the spin-independent parts of two 
Hamiltonians determining the particle and nucleus motion coincide formally, the obtained 
equations of motion are significantly different. The equation of particle and nucleus motion 
in the cylindrical system of coordinates includes the centrifugal force and also determines 
the effect of the mutual influence of particle motion in the directions e p and e^. This effect 
has a kinematic nature and, in particular, is manifested in the appearance of the oscillating 
force in the azimuthal direction during oscillatory motion in the radial direction and vice 
versa. Passing to the classical limit makes it possible to establish that the equations of mo¬ 
tion for the kinetic momentum and spin in the two considered systems of coordinates agree 
completely with each other. As an example demonstrating the correctness of the description 
of physical phenomena using the equations obtained in this paper, we have derived a well- 
known relation between the angular velocities of the spin and momentum vector rotations 
during planar channeling in bent crystals. 

Using methods developed in gravitation theory makes it possible not only to provide 
the correct description of particles and nuclei with the AMM and the EDM in cylindrical 
coordinates, but also to determine the physical nature of the main specific features of such 
a description. The analogy between the cylindrical system of coordinates and the rotating 
reference system demonstrated previously in [20] within the framework of the classical ap¬ 
proach appears paradoxical at first glance. Although the cylindrical system of coordinates 
and the rotating reference system belong to planar spacetime manifolds, the structures of 
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the metric tensor in these systems are significantly different. The main physical properties 
of the rotating reference system are determined by off-diagonal components of the metric 
tensor g 0i , and the metric is stationary and nonstatic. As opposed to this, the metric of 
the cylindrical system of coordinates (j9]) is static and has a single nontrivial component. 
However, the calculation of gravitoelectromagnetic fields explains the indicated analogy ex¬ 
haustively. In both cases, there is no gravitoelectric field, and the gravitomagnetic field & is 
uju° in the rotating reference system and (u^ / p)e z in the cylindrical system of coordinates. 
If the condition u ) = cue z = [u^ / (pu°)]e z is satisfied, these quantities are equal to each other, 
and the dynamics of particles and their spins in the two systems becomes identical. This 
condition can always be satisfied for an individual particle, but it cannot be satisfied for 
an ensemble (beam) of particles with different momenta. The possibility of describing an 
ensemble of particles with an arbitrary distribution over momenta is an important advantage 
of the cylindrical system of coordinates. 

Thus, in this paper, we presented the general quantum-mechanical description of rel¬ 
ativistic spin 1/2 particles and nuclei channeled in bent crystals by using the cylindrical 
system of coordinates. The results of our study can be used to solve concrete problems of 
channeling theory. 
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